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Abstract. The opacity of spiral galaxy disks, from counts of distant galaxies, is compared to HI column densiti es. The opacity 
measurements are calibrated using the "Synthetic Field Method" from Gonzalez et al. ( 1998); Holwerda et al. 1 2005a). When 
compared for individual disks, the HI column density and dust opacity do not seem to be correlated as HI and opacity follow 
different radial profiles. To improve statistics, an average radial opacity profile is compared to an average HI profile. Compared 
to dust-to-HI estimates from the literature, more extinction is found in this profile. This difference may be accounted for by an 
underestimate of the dust in earlier measurements due to their dependence on dust temperature. Since the SFM is insensitive to 
the dust temperature, the ratio between the SFM opacity and HI could very well be indicative of the true ratio. Earlier claims 
for a radially extended cold dust disk were based on sub-mm observations. A comparison between sub-mm observations and 
counts of distant galaxies is therefore desirable. We present the best current example of such a comparison, M51, for which the 
measurements seem to agree. However, this remains an area where improved counts of distant galaxies, sub-mm observations 
and our understanding of dust emissivity are needed. 



1. Introduction 

The relationship between gas and dust in spiral disks has been 
the focus of many observational studies often combined with 
efforts to characterize the chemical composition. Dust plays an 
important role in the energy and chemistry budgets of a disk, as 
is evident in our own Galaxy. The question comes simply down 
to whether the dust is distributed as the stars, which produce it, 
or like the gas, such as atomic hydrogen, which is dynamically 
coupled to it. Furthermore, some fraction of the gas which is 
associated with the dust will be in the form of cold dark molec- 
ular clouds, which may be of very high opacity. Cold dust can 
in principle be detected by either sub-millimeter emission or 
through the extinction of a background source. 

A first attempt to characterize cold du st clouds and their 
relation to HI in a nearby galaxy, was by iHodgel dl980t) : his 
figures 8 and 9 show the radial distribution of the number of 
dark clouds and HI column density along the minor and major 
axes of M3 1 . No correlation between the two tracers was found. 

In recent years, the Infrared Space Observatory (ISO), and 
the Sub-Millimeter Common User Bolometer Array (SCUBA) 
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on the James Clerk Maxwell Telescope have detected and 
mapped the cold dust component of spiral disks. This cold 
component has been found at mu c h larger radii in galaxy disks 
than th e warm dust. lAlton et ail dl998bl) and iTrewhella et alJ 
(2000) found evide nce from ISO for cold dust at larger radii. 
lAlton et all dl998al 1200 0b) studied the distribution of emis- 
sion in NGC 89 1 an d also found evidence for a cold dust disk. 
lAlton et all d2000al) found evidence for dust outside the opti- 
cal disk in NGC 660. In addition, iBianchi et ail d2000l) and 
lAlton et all d2004l) found a correlation between CO and 850 
micron emission in NGC 6946 and interpreted this as evi- 
dence for a correlation between molecular hydrogen and cold 
dust. The relation between gas and dust has been explored us- 
ing the SCUBA array and synthesis ma pping of the atomic 
hydrogen (HI). Re cent SCUBA results dStevens etaD 12005: 
iThomas e t alJ200 4|) are p resented in relation to HI column den- 
sity. [Stevens et al. (2005) gives the ratio between gas and dust 
masses and lThomas et alJ d2004l) compare the radial extent of 
dust and HI and find similar scalelenghts for both. 

Spectral ener gy distribution mod e ling (SEP) of edge- 
on galaxy disks dPopescu et al.l 2000: M isiriotis et al.l 12001 : 
Pop escu & Tuffsl2002l) indicate that the dust's emission is pow- 
ered by HII regions in the case of warm dust (A < 100 fjm.) and 
diffuse stellar radiation in the case of cold dust (A > 100 fim.). 
The infrared flux is dominated by dust illuminated by nearby 
stars and an cold dust is only visible in the sub-mm part of the 
spectrum. Therefore, a cold dust component is relatively more 
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prominent at higher radii where the dust clouds are not illumi- 
nated as much by the stellar population. 

Parallel to this observational effort to characterize the 
emission, there is an effort to characterize the extinc- 
tion in disks using known background sources. Two types 
of known background sources are in use: occulted galax- 
ies and the number of distant galaxies. The occulting 



galaxies retrieved (N) can be found and we fit the following 
equation to this: 



A = - 2 - 5CI "'[Zo 



(1) 



galaxy technique has been exhausted on the rare nearby 


pairs (White & Keel 


1992; 


Andredakis & van der Kruit 1992; 


Domingue et al. 


1999; 


White etai] 


2000; 


Domingue et al. 


2000; Keel & White 


2001 


albl). [Domineue et al.1 (Il999|) com- 



pared the extinction and dust emission in their pairs and found 
reasonable agreement between dust masses and no need for an 
extremely cold (T < 10K.) component. 

The present paper is one of a series using the number of dis- 
tant galaxies seen through the foreground disk in Hubble Space 
Telescope (HST) images as an extinction probe. The identified 
number of distant galaxies suffers from crowdi ng and confu- 
sion ef fects. To calibrate this observed number, Gon zalez et alJ 
Jl998h developed the "Synthetic Field Method" (SFM). A se- 
ries of synthetic fields is constructed, with the original science 
field to which a dimmed deep field is added. From the relation 
between detected added background galaxies and the dimming 
of the deep field, the average dimm ing of the science field can 
be inferred.lHolwerda et alJJ2005al) automated this method and 
IHolwerda et alJ J2005bl) reported on the radial opacity profiles 
of a sample of nearby galaxies. 

In this paper, t he radial opacity profiles from 
IHolwerda et ail j2C)05b) are compared to HI surface density 
profiles from the literature for the same galaxies. The inferred 
ratios and radial profiles are compared to those obtained from 
sub-mm observations and others in the literature. 

This paper is organized as follows: section [2 gives a brief 
description of the "Synthetic Field Method". The HI surface 
density radial profiles, radial opacity profiles and their ratios 
are presented in section [5] Radial opacity-to-HI profiles aver- 
aged for the whole sample are compared to values from the 
literature in section|3] Section[5]compares the sub-mm profiles 
from the literature to our opacity profiles. The conclusions are 
summarized in section |6] with a view to future work in section 

2. The "Synthetic Field Method" 

The number of distant galaxies seen through a spiral disk does 
not only depend on the level of extinction in the disk but also 
on the crowding and confusion by the objects in the fore- 
grou nd disk. To calibrate t he e ffects of crowding and confu- 
sion, jGonzile^^OC^l^ and lHolwerda et afl ( 2005a) devel- 
oped the "Synthetic Field Method" (SFM). This method con- 
sists of several steps. First the number of distant galaxies seen 
in the science field is identified. The selection is based on ob- 
ject characteristics and color and visually checked. Secondly, 
a series of synthetic fields is constructed. These are the origi- 
nal science field with a Hubble Deep Field added, dimmed to 
mimic dust extinction. Thirdly, the added distant galaxies are 
identified in these synthetic fields. A relationship between the 
dimming of the synthetic field (A) and the number of added 



where C characterizes the crowding and confusion for this par- 
ticular science field and No is the number of galaxies expected 
in the case of no dimming by dust. Substituting N by the actual 
number of galaxies found in the science field, equation^yields 
the average opacity for the field. The counts of distant galaxies 
are done in I band images and hence we report the opacities as 
A,. 

Cosmic variance in the number of distant galaxies in a 
given field adds an extra uncertainty to the number of distant 
galaxies found in the science field. As a result, opacity mea- 
surements in a single WFPC2 field or section thereof have high 

uncertainties associated with them. To combat this, we have 

■ ML 
automated this method ( Holwerda et al. 2005a) and applied it 

to a si zeable sample of archival WFPC2 fields JHolwerda et alJ 
2005b). In this series of papers, we have explored the relations 
of disk opacity with radius (Holwerd a et alJl200 5b). the sur- 
face brightness (IHolwerda et alJl2005d) and HI in this paper. 
Holwerd a et alJ (|2005dl) explore the li mitations of this method 
as predicted by Gonzale z et alJ (12003), concluding that the op- 
timal foreground disk distance is somewhere between 5 and 30 
Mpc. 

3. Individual galaxies: radial HI and opacity 
profiles 

The number of distant background galaxies as a function of HI 
column density is best directly measured using an overlay of 
the HI column density map on an HST field. See for example 
the analysis by Cuillan dre et alJ (1200 II) of a ground-based field 
in M31. However, since the HI column density maps are not 
easily available and not uniform, the relation between radial 
profi les of HI and opacity is used. 

IHolwerda et ail |2005b) present radial profiles for individ- 
ual galaxies and composites of fields (their Table 3). HI surface 
density profiles were taken from the literature for the subset of 
our sample for which these were available (Table^. HI surface 
density profiles of the galaxies we re extracted from the litera- 
ture using the DEXTER program ( Demleitner et al. 2001) and 
rescaled to express radius in ^25 JdT^uToTrteurset al Jl 1 99 ll) 
and surface density in M pc~ 2 . 

In Figure ^3 and ^] the individual radial surface density 
profiles of HI and the opacity profiles are shown. The opacity 
was determined for the sections of the WFPC2 fields corre- 
sponding to radial intervals of 0.25 /?25- Both the opacity pro- 
files and the HI surface density profiles of spiral disks display 
a variety of shapes. However, in general, the HI profile peaks 
somewhere in the disk and flattens out or dips near the galaxy's 
center. The opacity profiles show a gradual rise towards the 
disk's center. It should be noted that the HI profiles are the az- 
imuthally average for the entire disk, while the opacity profile 
is derived from the smaller section of the disk corresponding to 
the HST image i.e. these profiles are for the same spiral galaxy 
but not determined from the same section of the disk. 
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To see if there is correlation between HI surface density and 
disk opacity, the values of the profiles in Figures ITal and II blare 
plotted in Figure|2 Averaged over radial intervals of 0.25 R25, 
there seems to be no correlation between HI column density 
and dust opacity in a spiral disk. The lack of a relation may be 
explained by the fact t hat the opacity p rofiles generally rise in 
the center of galaxies dHolwerda et al.l2005bl) . while radial HI 
surface density profiles often show a drop in the galaxy's center 
(See also Fi gures ITal and ITbl . Alternatively, part of the hydro- 
gen gas may be associated in molecular clouds, undetected in 
HI observations. 



4. Dust-to-HI ratio 

To reduce the uncertainties in the opacity measurement, the 
galaxy counts from several fields must be combined. This can 
be done e.g. pe r Hubble type or for the entire sample (See 
IHolwerda et al] (E"0Q5b) for this type of analysis). 

In Figure [5] the average Hl-to-opacity plot based on our 
entire sample (Table [0 is compared to Hl-to-opacity measure- 
ments from the literature. To obtain this Hl-to-dust profile, the 
opacities were re-derived and the HI profiles averaged. First, 
all the counts of distant galaxies, from all the synthetic and 
science fields, were combined. Then the opacity was derived 
again from the combined counts in radial intervals, using equa- 
tion n The HI surface density profiles from figures [Talandllbl 
were combined by adding mass and surface contributions of 
each profile in a radial section and taking the ratio of the sums. 

An average opacity taken over many galaxies does not suf- 
fer from poor statistics as the individual profiles do, but small 
scale variations are smoothed out. The radial profile of dust- 
to-HI remains relatively constant for most of the optical disk 
and turns upward beyond that. At this point, the HI profiles di- 
minish but the opacity profiles remain more or less constant. 
The radial profile becomes more uncertain beyond R25 as the 
WFPC2 fields used for the SFM analysis were centered on the 
optic al disk (See fo r more o n the selection of our original sam- 
ple Holwerda et al . (2005b) .). For illustration, only two galax- 
ies have enough counts for a reasonable SFM measurement be- 
yond the R25 ■ The trend of the dust-to-HI ratio beyond the R25 
is therefore correspondingly uncertain. 

The comparison data from the literature are b ased on a 
series of techniques to determine the dust content. Ilssa et all 
(1990) compile dust-to-HI measurements from the litera- 
ture for a few n earby galaxie s and normalize these to the 
Galactic value of Boh lin et al.l dl978l) and a distan ce of 0.7 
R^ fro m their galaxies's center. These values from Ilssa et alJ 
( 1990) a re derived from the line ratio between Ha and HyS. 
Mav va & Reng arai a infer average dust-to-HI profiles 

from IRAS emission at 60 micron which traces only the warm 
dust. They give a range of conv ersion values for th e optical 
depth at 60 micron and in V. iBoissier et alJ d2004l) present 
radial dust-to-HI profi les for 6 galax ies deduced from the 
FIR/UV ratio. Cuillandr eetail (2001) derive a relation be- 
tween HI column density and stellar reddening for M3 1 . The 
average Hl-to-dust ratio from the sub-mm observations of disks 
by IStevens et al.l (120051) is also indicated. Their total sub-mm 



flux den sity was converted using the ratio found bv lAlton et alJ 
( 1998b) and an estimate of the disk area using the D25. 

The comparison to other results in Figure|3]raises the ques- 
tion why the opacity measurement from the number of distant 
galaxies results in a dust-to-HI ratio that is an orde r of mag- 
nitude higher than the nearest recent estimate dBoissier et alJ 
2004). There are several effects which could have affected the 
average dust-to-HI profile from counts of galaxies. (1) The 
opacity measure from the counts of galaxies is averaged over a 
series of Hubble types and HI profiles. (2) The older HI surface 
density profiles may underestimate the HI column density, es- 
pecially the larger scale structure due to lack of short-spacing 
information. (3) The opacities are predominantly determined 
for fields with a spiral arm while the HI profiles are averages 
over the whole of disk. This may have lowered the HI column 
density profile with respect to the measured opacity. (4) The av- 
erage opacity measurement from counts of galaxies indicates 
the covering factor of dark clouds. One could argue that this 
covering factor is an overestimate as the distant galaxy does 
not need to be completely covered in order to be dropped from 
the counts. All these effects could have lowered the dust-to-HI 
profile presented here. 

The other profiles also suffer from various different system- 
atics. The dust-to-HI ratio estimates from the literature rely on 
the light of the disk itself to estimate the dust content. This in- 
herently biases the measurement to lower extinction values, as 
not the entire height of the dust disk is responsible for the ob- 
served extinction and the measured light is biased towards low- 
extinction lines of sight. The different manner in which the dust 
content was es timated may also play a role. The values from 
Ilssa et all dl990l) are for very specific parts of the disk, e.g. the 
HII regions, and these were scaled to a single radius. Still, the 
ratio is similar to those found by other auth ors for the whole of 
the disk. The infrared measures of Mavva & Rengarai an ( 1997) 
(A = 60 /urn.) show substantially more HI for a given amount of 
dust. Most likely this is the effect of their observational tech- 
nique which is more sensitive to the warmer component of the 
dust in the spiral disk, which dominates the infrared emission. 
The fact that infrared emission onl y detects 10 to 20 % of the 
dust in a disk was also found by iDevereux & Yound (|l990) 
from 60 and 100 ymi. emission obs erved with IRAS. By tak- 
ing the ratio of UV and far-IR flux, IBoissier et al.l d2004l) can 
characteriz e better the total dust content of a disk. The average 
value from St evens et alJ d2005l) is again closer to this paper's 
profile. 

Th e progression from the values of |Mavva & Rengarajanl 
19971) to those of IBoissier et alJ (120041) and IStevens et alJ 
2005) indicates that the more sensitive the dust indicators be- 
come for the colder dust component in a disk, the more dust is 
detected in r elation to the HI. It has also been found by sev- 
eral authors dAlton et alJll998bt IPooescu & Tuffs! Eo02) that 
the bulk of the dust reveals itself when more sensitive dust in- 
dicators are used. Taken in this context, the SFM estimate (the 
thick line in Figure |5J of dust content of disks, which is inde- 
pendent of dust temperature, may reflect the actual total dust- 
to-HI ratio. 



4 



Holwerda et al.: dust emission, extinction and HI 



5. Comparison to SCUBA profiles. 

The SCUBA detector on the JCMT detects sub-mm emission 
from cold dust in spiral disks. In order to verify whether or not 
this cold dust is responsible for the opacity measured with the 
SFM, a direct comparison should be made. Four galaxies in our 
s am ple were mapped in t he 850 uan band by di f ferent authors 
jAltonetalJl200n 120021 iMeiierink etaD 120051: IStevens et ail 
120051) . However, the individual opacity measurements and the 
conversion from 850 micron flux to an optical depth do not al- 
low for a good comparison. The exact emissivity o f dust g rains 
at these wavelengths may still be underestimated jAlton et alJ 
2000jJ l2004l lDasvraetai1l2005l) . To illustrate a comparison 
between a SCUBA map and our counts, we p resent the best 
curren t example: the exponential disk found by Meii erink et alJ 
J2005I) for M5 1 and our distant galaxy counts from two WFPC2 
fields. 



5.1. M51 

IMeiierink et al.l J2005I) present a detailed map of the 850 mi- 
cron emission, as well as an exponential optical depth profile 
for the d ust disk. Figure Pflshow s the 850 micron emission pro- 
file from Meiierin k et alJ (120051) translated to an opacity value 
in V and I and the opacity measurements of the two WFPC2 
fields and their average. Our WFPC2 fields predo minant ly miss 
the sections of the disk for which Meiierink et al. (2005) report 
additional flux from the spiral arms. However, a proper compar- 
ison can be done when the sub-mm observations and the counts 
of distant galaxies are from a similar sized field. For now we 
can only conclude that the points and the profile seem to agree. 



5.2. Comparing to sub-mm 

A good comparison between opacity and sub-mm flux can be 
obtained with more and better sub-mm maps and sufficient 
counts of distant galaxies. The new SCUBA-2 instrument on 
the JCMT promises to facilitate this mapping of nearby galax- 
ies. However the expected emissivity of dust grains must then 
also be known for an accurate comparison. An advantage of the 
counts of distant galaxies is that they can be extended to larger 
radii and lower opacity values and can independently verify the 
profile found from sub-mm observations. 

Both observational techniques can be used to explore the 
radial exte nt of dust in spiral di sks. Our previous compos- 
ite profile faolwerda et alJl2005bl) already pointed to and ex- 
tend ed dust disk as well as several sub-mm observ ations 
llNelson et al ill 998t lAlton et alJIl 998ct iTrewhella et alJkoOO: 
Pope scu & Tuff si 12003). Currently the extent is taken to be 
somewhere between the optical and the HI scale. A more ac- 
curate determination of that can be achieved with either more 
sensitive sub-mm observations or future galaxy counts or both 
combined. 



6. Conclusions 

In summary, we draw the following conclusions from the com- 
parison between opacity from the number of distant galaxies 
with HI profiles: 

1 . The HI and opacity profiles of individual fields do not seem 
to correlate well. This may be due to the rise of opacity in 
the center of a disk while there is often a dip in the atomic 
hydrogen surface density profile there. However, high un- 
certainties plague the single field opacity measurements 
(Figures[liITbland|3. 

2. The dust-to-HI ratio depends on the tracer used for the dust 
content. When a tracer sensitive to colder dust is used, a 
higher dust-to-HI value is found (Figure[3}. The SFM opac- 
ity, which is independent of the dust temperature, could 
very well point to the true dust to HI ratio in spiral disks. 

3. A direct comparison between sub-mm emission and SFM 
opacity is problematic as few galaxies in our sample have 
been observed at these wavelengths. It is possible that both 
methods trace the same component (Figure^ but more data 
and a better understanding of dust emissivity is needed for 
confirmation. 

7. Future Work 

The relation between sub-mm emission and opacity is an im- 
portant area for future work. A more accurate comparison be- 
tween the SFM opacity profile and the SCUBA emission pro- 
file of M5 1 should become possible in the near future with the 
HST/ACS observations by the Hubble Heritage Team. A simi- 
lar comparison could be made for M101, provided such a large 
solid angle could be successfully scanned by the new SCUBA- 
2 or its successors. 

By comparing the SFM opacity to HI column density aver- 
aged over a radial interval, small variations are smoothed out. 
A more direct approach is to compare SFM opacity directly 
in contours of HI column density. This requires uniform HST 
imaging of a nearby galaxy and a deep and uniform HI map as 
well. Such data exists for M101 and M51 and the relation be- 
tween HI and dust could then be characterized more accurately 
and out to a larger radius. 
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Table 1. 



Galaxy 


Hubble Type 


Instrument 


^25 

arcmin. 


Reference 


NGC 925 


SAbd 


WSRT 


10.47 


Weversetal. (1986) 






VLA 




Pisanoetal. (19981 


NGC 1365 


SBb 


VLA 


11.22 


Jorsater & van Moorsel (1995) 
- 


NGC 1425 


SBb 


- 


5.75 


NGC 1637 


SAB(rs)c 


- 


3.98 


- 


NGC 2541 


SAcd 


WSRT 


6.31 


Broeils & van Woerden (1994) 


NGC 2841 


SAb 


WSRT 


8.13 


Bosma (19811 


NGC 3198 


SBc 


WSRT 


8.51 


Wevers et al. (19861 


NGC 3319 


SB(rs)cd 


WSRT 


6.17 


Broeils & van Woerden (1994) 


NGC 3351 


SBb 


- 


7.41 


- 


NGC 3621 


SAc 


- 


12.3 


- 


NGC 3627 


SAB(s)b 


- 


9.12 


- 


NGC 4321 


SABbc 


WSRT 


7.41 


Warmels (1988) 






VLA 




Cavatte et al. (1994) 


NGC 4414 


SAc 


VLA 


3.63 


Thornlev & Mundv (1997) 


NGC 4496A 


SBm 


- 


3.98 


- 


NGC 4527 


SAB(s)bc 


- 


6.17 


- 


NGC 4535 


SABc 


WSRT 


7.08 


Warmels (1988) 






VLA 




Cavatte et al. (1994) 


NGC 4536 


SAB(rs)bc 


- 


7.59 


- 


NGC 4548 M91 


SBb 


WSRT 


5.37 


Warmels (1988) 






VLA 




Cavatte et al. (1994) 


NGC 4559 


SAB(rs)cd 


WSRT 


10.72 


Broeils & van Woerden (1994) 


NGC 4571 


SA(r)d 


WSRT 


3.63 


Warmels (1988) 








j.jy 




NGC 4639 


SABbc 


WSRT 


2.75 


Warmels (1988) 


NGC 4725 


SABab 


WSRT 


10.72 


Weversetal. (1986) 


NGC 6946 


SAB(rs)cd 


VLA-D 


11.48 


Tacconi& Youna(1986) 


NGC 7331 


SAb 


WSRT 


10.47 


Bosma (1981) 


M51 


SA(s)bc 


WSRT 


5.61 


Bosma (1981) 


M81 


SA(s)ab 


WSRT 


13.46 


Cavatte et al. (1994) 
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Fig. la The HI and opacity profiles of individual galaxies in Fig. 2 The average opacity (A/) and surface density (Ehi) for 
our sample. A/ is expressed in magnitudes, the atomic hydro- radial annuli of 0.25 R25 in individual galaxies in our sample 
gen surface density (Shi) in solar masses per square parsec (See Figures [la]and[Tb). 
(M pc~ 2 ). The HI surface density profiles often do not have 
uncertainties reported with them. See Table ^for the original 
reference for the HI surface density profiles. The most recent 
profiles are used in all cases. NGC 4414 has two WFPC2 fields 
associated with it. 
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Fig. 3 The radial dust-to-H I profiles from Issa et alJ {l990), 
iMavva & Rengaraian|(ll997h . lBoissier et alJd2004 and this pa- 
per (thick l ine, uncertainties are dashed). The symbols are the 
values from llssa et alJ l ll990m Ha7HS ratio). Not e how these are 
higher for the same galaxy as the valu es from Boissier et al. 
(2004). The three bottom lines are the Mavva & Rengaraian 
( 1997) (IRAS 60 fim emission) average ratio, for the different 
conversion factors between the 60 fim optical depth and the V 
band optical depth (j v / T sO). The curves in t he middle are 
the dust-to-HI profile from Boissier et al. (2004)(UV/FIR ra- 
tio), for their 6 gala xies, converted to I band extinctions using 
Bosel li et al. (2003). The arrow right is the average ratio from 
Stev ens et alJ 1120051) (850 ymi emission), the combined warm 



Fig. 4 The profile of M5 1 fro mlMeiiermk et al l d2005l) (dashed 
lines) and the opacities from lHolwerda et alJ J2005b) (points). 
Solid points are the individual WFPC2 fields, the open triangles 
are the derived from combined WFPC2 fields. The errorbars 
shown are for the combined opacity and the radial bin size is 
0.25 /?25- The top axis shows the radius in R^S). Conversion 
from the SCUBA profile to an opacity in V used two times the 
Galactic emissivity at 850 micron and the Galactic reddening 
law was used to convert the opacity profile in V to one in I. 



and cold dust over HI. 



